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The radial collapse of a homogeneous disk of collisionless particles can be solved analytically in Newtonian gravitation. To 
solve the problem in general relativity, however, requires the full machinery of numerical relativity. The collapse of a disk 
is the simplest problem that exhibits the two most significant and challenging features of strong-field gravitation: black hole 
formation and gravitational wave generation. We carry out dynamical calculations of several different relativistic disk systems. 
We explore the growth of ring instabilities in equilibrium disks, and how they are suppressed by sufficient velocity dispersion. 
We calculate wave forms from oscillating disks, and from disks that undergo gravitational collapse to black holes. Studies 
of disk collapse to black holes should also be useful for developing new techniques for numerical relativity, such as apparent 
horizon boundary conditions for black hole spacetimes. 
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I. INTRODUCTION 

The simplest example of gravitational collapse is that 
of a homogeneous sphere of particles initially at rest. 
This collapse solution is analytic both in Newtonian grav- 
ity and general relativity. In general relativity, this so- 
lution is known as Oppenheimer-Snyder collapse (the 
solution is a "piece" of a closed Friedmann universe) . Be- 
cause of Birkhoff's theorem we know that this solution is 
nonradiating. Both the particle motion and the gravita- 
tional field are radially symmetric, i.e., functions of one 
spatial variable. 

The radiating problem which is the simplest analogue 
to Oppenheimer-Snyder collapse is that of an axisymmet- 
ric, infinitely thin disk of particles initially at rest. This 
case is simple because the particle motion still depends 
on only one spatial coordinate, although the gravitational 
field now depends on two. If the disk is constructed by 
squashing a homogeneous sphere into a pancake, keep- 
ing the density homogeneous, then the solution is still 
analytic in Newtonian theory 

However, to solve the problem for relativistic gravita- 
tion requires the full machinery of numerical relativity, 
and has not been addressed until now. Indeed, the dy- 
namical properties of a disk of collisionless particles has 
never been studied in general relativity, although it has 
been extensively treated in Newtonian theory || . In this 
paper we tackle the collapse of an axisymmetric collision- 
less disk of particles in general relativity theory. Particles 
in an axisymmetric disk can also have angular motion, 
but because of conservation of angular momentum this 
motion is not dynamical. We consider cases in which the 
disk particles are initially at rest and also in which they 
have initial angular motion. In the latter case we focus 
on disks with total J = 0, i.e. with equal numbers of co- 
and counter- rotating particles. 

In Newtonian theory, it is known that equilibrium disks 



supported against collapse by rotation alone are unsta- 
ble to ring formation |5|. The disk can be stabilized by 
"heating" the disk, that is, converting some of the or- 
dered rotational energy into random "thermal" motion 
. We explore here whether or not a similar result holds 
in general relativity. 

Our study of disks is geared to analyze two relativis- 
tic phenomena that do not arise in Newtonian theory: 
collapse to black holes and the generation of gravita- 
tional waves. Since this is the simplest wave generation 
problem, disk collapse provides a useful proving ground 
for testing codes designed to treat gravitational radia- 
tion in general relativity. Once we can evolve disk sys- 
tems accurately in general relativity, we should be able 
to compute the gravitational field from any axisymmetric 
source, since the equations for the field are essentially the 
same as those for general axisymmetric sources. The ma- 
jor problems associated with general axisymmetric col- 
lapse are all contained in this test case: formation and 
evolution of black holes and propagation of gravitational 
waves. 

We have previously studied disk collapse in the context 
of nonlinear scalar gravitation, where many of the tech- 
niques for handling infinitely thin disks were developed 
Q . The basic code for evolving axisymmetric spacetimes 
in general relativity has been discussed in a number of 
papers [BlHJ^l ■ We will refer extensively to the equations 
in Ref. Bin the discussion that follows. The organiza- 
tion of this paper is as follows. In Secjn] we present the 
equations for an evolving, general relativistic disk. In 
Sec. Ill we discuss analytic test problems including New- 



tonian and relativistic disk solutions. In Sec. IV we give 
the results of our numerical calculations. 
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II. BASIC EQUATIONS 

The particles comprising the disk are assumed to in- 
teract exclusively by gravitation, i.e., they obey the rel- 
ativistic collisionless Boltzmann equation (Vlasov equa- 
tion) . Accordingly, we have constructed a numerical code 
that solves Einstein's equations for the gravitational field 
coupled to matter sources obeying the Vlasov equation. 
This is the same mean-field, particle simulation code de- 
scribed in Refs. |]|| to study nonspherical gravitational 
collapse. The code is designed to handle axisymmctric 
systems with no net angular momentum. The present 
version assumes equatorial symmetry. We solve the field 
equations in 3 + 1 form following Arnowitt, Deser, and 
Misner ||. We use maximal time slicing and quasi- 
isotropic spatial gauge in axisymmetry. The metric, writ- 
ten in spherical-polar coordinates is 



ds 2 



l dt 2 + A 2 (dr + I3 r dty + A 2 r 2 {d9 + (3"dt) 



B 2 r 2 sin 2 9 dcp 2 . (1) 



The matter satisfies the relativistic Vlasov equation, 
which we solve by particle simulation in the mean grav- 
itational field. The basic code is identical to the one 
described in Refs. ||. The key equations and definitions 
of variables are given in the Appendix of Ref. || . (These 
were generalized to include net rotation in Ref. &.) 



A. Jump conditions 

The disk matter source affects the metric via jump con- 
ditions in the field equations across the equatorial (disk) 
plane. These jump conditions replace the usual matter 
source terms that appear in the field equations. A deriva- 
tion and numerical implementation of such a jump condi- 
tion was presented in Ref. for the simple case of scalar 
gravity. There the governing equation is a nonlinear wave 
equation with a matter source on the right-hand-side. 

As an example of how the jump conditions may be 
derived in 3 + 1 general relativity, consider the two mo- 
mentum constraint equations (A4) and (A5) in Ref. 0: 
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Since the particles are confined to the equatorial plane 
9 = 7r/2 where (3 e — 0, the particle 4-velocity component 
satisfies u e = ug = 0. Hence So = 0. Equation (|^) 
integrated across the equator yields 



-dg( S m6K r 9 ). 
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where ± denotes 9 — it/2 ± e, e — * 0. Functions that 
are symmetric across the equatorial plane, such as k r r 
and K^, are continuous there. Hence Eq. (|4|) reduces to 
= 0. Now consider Eq. (|J). Integrating it gives 



= 



S r r sin I 



- -k 
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Since k@ is antisymmetric across the equator, Eq. || gives 



k r n \ + = -ki\- = - 



S r rsin9d9. 



(6) 



Similarly, integration of the Hamiltonian constraint equa- 
tion (Eq. (A6) of Ref. §) leads to 

1 1 1 f + 

-sinOipa] + = -—ipr},B ~ rrr / P*rsin9d9. (7) 
r At sip J_ 

Integrating the lapse equation (Ref. || Eq. (A7)) gives 



- sm9(aip) g\ + = ——aipr) , 
r 4r 



\^ [ (p* + 2S)r sin 6d9. 
8 B J _ 



(8) 



The boundary condition Eq. ^ is used to set the value of 
Kg all along the equatorial plane. In the vacuum, outside 
of the equatorial plane, Kg is determined by integrating 
the evolution equation, Eq. (A3) of Ref. ||, as usual. 

When finite differencing the Hamiltonian constraint 
(Eq. A6 of Ref. |(|), the derivative terms ip o and rj t g ap- 
pear in exactly the combination as in Eq. (Jt]). The only 
place where the matter source term p* appears in the 
Hamiltonian constraint is through this boundary condi- 
tion. Equation ^ is used in an analogous fashion for the 
lapse equation (Ref. |(| Eq. A7). 

The dynamical equation for r\ (Ref. M Eq. A2) and 
the shift equations (Ref. § Eqs. A8-A9) for (3 r and 
remain unchanged. Note that 77, (3 r and (3^ are metric 
coefficients and thus must be continuous across the equa- 
tor. 



B. Matter sources 

The geodesic equations of motion for the collisionless 
matter particles are given by Eq. (A10-A16) of Ref. § 
with the following simplifications: ug — and 9 = tt/2. 
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Hence, as is spherical symmetry, only the radial motion 
is dynamical for an infinitely thin disk. 

The particles are binned in annuli to determine the 
source terms for the field equations. Equations (A17- 
A21) of Ref. § lead to the disk sources f§: 



a = I p r sin 
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Here m is the particle rest-mass related to the total rest- 
mass Mo by m = Mo/N with N the total particle num- 
ber. We obtain Mq from 



Mo = / ao2irrdr, 



where 



co = / Par sin ( 
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and where pg is the rest-mass density. 



C. Hydrodynamical disks 

In the special case in which concentric shells of parti- 
cles do not cross, collisionless matter may be treated as 
hydrodynamical dust. Thus, as an alternative to inte- 
grating geodesic equations followed by binning of parti- 
cles, one could integrate the equations of relativistic hy- 
drodynamics. The hydrodynamical equations have the 
disadvantage that they are PDE's and not ODE's like 
the geodesic equations. However, they have the advan- 
tage that they produce intrinsically smooth source pro- 
files, unlike particle descriptions which are stochastic. 

The basic equations of relativistic hydrodynamics in 
3+1 ADM form are given in, e.g., Ref. [^0|. For a cold 
axisymmetric disk they reduce to the continuity and ra- 
dial Euler equations: 

<9 t er + -d r (rcr v r ) = 0, 
r 

9 t S r + -d r (rJ: r v r ) = -ad r a + Z r d r (3 r + Ed r In Aa, (15) 
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III. ANALYTIC SOLUTIONS AND TESTS 

Before we consider numerical solutions of the dynam- 
ical equations for disks and their gravitational fields, we 
give here some analytic results that will serve as code 
checks and initial data for our evolutions. 



A. Oscillating Newtonian disks 

As discussed in Ref. 0, there exists a complete an- 
alytic solution that furnishes a good test of a numeri- 
cal disk code in the weak field, slow- motion limit. The 
solution describes an oscillating homogeneous spheroid 
in Newtonian gravitation in the disk limit (eccentricity 
e — > 1). The surface density of such a disk of mass M 
and radius R is 
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Start with the equation of motion for the semi-major 
axis R of an oblate homogeneous spheroid (e.g. Eq. 5.6 
of Ref. |l|). Take the limit e -> 1 and find 



R = - 



3tt GM 



R 3 ' 



(20) 



where h is the conserved angular momentum per unit 
mass of a particle at the surface. Since the motion is 
homologous, the radius of each particle satisfies a similar 
equation. Choose h to be a fraction £ of the equilibrium 
angular momentum h = (3nGMRo/<i) 1 / 2 . Set 



R = R X(t). 



(21) 



Then the radius r of each particle satisfies 

r - r X(t), (22) 

where ro is the initial radius. Substituting Eq. ( ^l|) into 
Eq. (20) we see that X satisfies the familiar equation 
of an elliptic orbit for a particle with specific angular 
momentum h c s — ^(M/i?^) 1 / 2 around a fixed central 
mass M e g — Z-kM/A. The parametric solution for X(t) 
is 



(23) 
(24) 



X = a(l — e cos u), 
t = —{u-esmu) - —. 



where we assume X = 1 and X = at t = 0. In Eqs. Q25 
and (|24|), the semi-major axis, eccentricity and period are 
given by 
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The radial and tangential particle velocities are given by 
X 
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It is simple to derive the gravitational wave amplitude 
for an oscillating disk in the quadrupole approximation. 
In axisymmetry, in the absence of rotation, there is only 
one polarization and its amplitude is given by 



r/i-i 



:Izz sin 



(28) 



where 



!zz = ~\ Jr 2 pd 3 x = ~ Jr 3 adr = -^MR 2 (29) 

Here we have used Eq. (|l9|). Using Eqs.(pp|) and ( pl| ) we 
find 
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B. Kalnajs Disk 

When £ — 1, the above Newtonian disk solution cor- 
responds to a uniformly rotating disk in dynamical equi- 
librium. As mentioned in the Introduction, such a disk is 
unstable to the formation of rings but can be stabilized 
by heating. Kalnajs Q has given an analytic prescription 
for constructing hot homogeneous disks in equilibrium. 
They all have the same surface density a and gravita- 
tional potential <!> as the cold disk in subsection A, but 
differ in the amount of random motion. In these models, 
the particles have an isotropic velocity distribution in a 
rotating frame that moves with angular velocity Q: 



v cos x = v 4> ~ 
v sin x — Vr ■ 



(31) 
(32) 



Here v is the magnitude of the isotropic velocity in the 
rotating frame and x is a random angle about the particle 
position in the disk plane. The distribution of particle 
velocities is given by 



f(v)vdv = 2-kK [< ax - v 



-1/2 



vdv. 



(33) 



Here 



{nl-n 2 ){Rl-r 2 \ 

hp 



(34) 
(35) 



and K is a normalization constant. Models in this family 
are parametrized by the ratio Q/Qq. Cold disks have 
£1 = f^o- Hot disks with O / Oo < 0.816 are stable against 
ring formation. 



C. Relativistic disk 



Here we construct a relativistic generalization of the 



cold homogeneous Newtonian disk described in Sec. Ill A 



This is useful for studying the dynamical behavior of 
disks in the strong field region, where almost nothing 
is known. 

Consider a disk of particles in circular equilibrium. 
The Hamiltonian equation (Ref. & Eq.(A6)) reduces to 



V 2 i/> = -2tt^-, 



(36) 



where we have temporarily restored a factor of 8-7T to the 
right-hand-side. We can obtain an analytic solution if 
we first cast Eq. (|36] ) into Poisson's equation by setting 
p*/ip = 2pw, obtaining 



V V = -47rp 



(37) 



and then equating pjv to the homogeneous Newtonian 
density profile for an oblate spheroid in the pancake limit 
El. We then find 



if) = 1 - $ 



N ■ 



(38) 



where $at is the Newtonian potential for a flattened 
spheroid El. In the disk interior, this yields 



3nM Nn lr 2 



4i? 



(interior). (39) 



The total mass of the disk M is related to the Newtonian 
mass M/v appearing in $jy by 



M 



Comparing Eqs. (|36|) and (|3 
density Eq. (g) is given by 

a = 2ipa N , 



2M N . (40) 
we find that the surface 

(41) 



where cat is the corresponding Newtonian density given 
by Eq. ©. 

The angular motion of a particle in an equilibrium disk 
is determined by setting du r /dt = u r = and using 
Eqs. (A.8) and (A.13) of Ref. f§. The result is 



(a, r /a)B 2 r 2 



B tr /B - a. r /a + l/r 



(equilibrium), (42) 
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where B = ip 2 is given by Eqs. (|3^) and (p9[). The lapse 
a is found from the maximal slicing [K\ = <9t-?Q = 0) 
condition, Eq. (A7) of Ref. |(| specialized to equilibrium, 
for which 77 = = KJj. The source term S appearing in 
the jump condition for the lapse equation is given by 



B 2 r 2 + ui 



(43) 



where we have used Eq. (A16) of Ref. ||. Because a 
and 11^ are interdependent, it is necessary to iterate the 
lapse equation and Eq. (f42|) . As an initial guess we use 

-1 , ^(exterior) 

The rest- mass surface density <7o defined in Eq.(|13j) is 
given by 



^0 



{l + ul/B 2 r 2 )V 2 



(44) 



from which the total rest mass can be computed using 
Eq. (gj). 

We can obtain a solution to the initial value equations 
for a nonequilibrium disk at a moment of time symmetry 
when the particles are moving at a fraction £ of their 
equilibrium velocity: 



u <t> = 



a) 



(45) 



The other equations remain unchanged. In the Newto- 
nian limit, this solution goes over to the cold oscillating 
disk of Sec. ETTA. 



An interesting nonequilibrium solution is the one in 
which £ = 0. This corresponds to the collapse of a homo- 
geneous disk in which all the particles are at rest initially. 
This is the pancake analogue of Oppenheimer-Snyder col- 
lapse of a homogeneous sphere. In the disk case, however, 
the solution is radiative and is not known analytically. In 
this case the total mass and rest mass are related by 



M = 



(see Ref. @). 



2M 



l + (l + 6^M /5i? ) 1 /2' 



(46) 



IV. NUMERICAL RESULTS 

In this section, we give examples of evolutions com- 
puted with our fully relativistic particle disk plus gravity 
code. When possible, we show gravitational wave forms 
and also make contact with analytic results. We also 
monitor quasi-local mass indicators; for the runs shown 
here the total mass is numerically preserved to within 
about 3%. In Table I we summarize the cases discussed 
in the text. 



FIG. 1. Gravitational wave form from a cold oscillating 
disk with an analytic Newtonian matter source. The numeri- 
cally generated and propagated gravitational wave form (solid 
line), extracted at a radius of AIM, is compared with the an- 
alytic prediction, Eq. (30) (dotted line). The disk has an equi- 
librium radius Ro/M = 30 and a velocity cutdown £ = 0.90. 
The wave amplitude h+ is dimensionless while r and t are in 
units of M. 



A. Oscillating Cold Newtonian disk 

As a check on our code, one would first like to simulate 
the coll apse o f the oscillating homogeneous disk described 
in Sec. Ill A . Unfortunately, such a cold configuration is 
unstable to ring formation. This is the familiar result for 
cold equilibrium disks with fl = flo and £ = 1 discussed 
above; we find numerically that it also holds for cold os- 
cillating disks. We can still employ this analytic model to 
check the field solver in our code provided we supply the 
unperturbed source function a analytically as a function 
of time. We then allow the code to solve the field equa- 
tions and thereby determine the gravitational radiation 
numerically from the metric and extrinsic curvature. 

Here we give results from a typical evolution of a ho- 
mogeneous disk with Rq/Mq — 30 and velocity cutdown 
factor £ = 0.9. This choice of radius is sufficiently large 
that the system remains essentially Newtonian through- 
out its evolution. We place the outer boundary at 
fmax/Mo — 500 and use a mesh with 30 radial zones in- 
side the matter, 220 outside, and 16 angular zones in the 
upper hemisphere. In Fig. [j]we compare the gravitational 
wave form extracted at an arbitrary fixed exterior radius 
of 42M using a gauge- invariant technique Jf| , with the 
analytic prediction Eq. (|30|). After an initial transient, 
the wave forms coincide very closely. The slight noise in 
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FIG. 2. Snapshots of the particle positions for the evo- 
lution of an oscillating Kalnajs disk. The initial radius is 
Ro/M — 20.0, velocity cutdown £ = 0.9, and velocity disper- 
sion parameter fi/fip = 0.8. Particle coordinates and time 
are in units of M. 



FIG. 3. Gravitational wave form from the oscillating 
Kalnajs disk shown in Fig. 2. The quadrupole wave form 
extracted at a radius of r/M = 21 (solid line) is compared 
with the quadrupole formula result, Eq. (28) (dashed line). 



the extracted wave form arises from interpolation error 
when the analytic source is mapped onto the numerical 
grid. 

The excellent agreement between the analytic and nu- 
merical results here gives confidence in the reliability of 
the numerical implementation of the jump conditions and 
the solution of the field equations. 



B. Oscillating Kalnajs disk 

To test the particle simulation aspects of our code, we 
set up stable equilibrium Kalnajs disks as described in 
Our code successfully holds these in equilib- 



Eq. 



with 



Sec. Ill B 



rium for several rotation periods. 

We next study the production of gravitational radia- 
tion from the oscillation of a hot nonequilibrium disk. 
We set up a Kalnajs disk with Rq/Mq = 20 and veloc- 
ity dispersion fi/fig = 0.8. We then induce collapse by 
reducing all particle velocity components by a cutdown 
factor £ = 0.9. This run employed 200 radial by 32 angu- 
lar zones and 12000 particles. Frames from the evolution 
are shown in Fig. ||. 

Although some of the collisionless matter from the hot 
disk escapes, most of the particles stay together and os- 
cillate homologously around the equilibrium radius. In 
Fig. ^ we compare the extracted quadrupole wave form 
with h + computed numerically from the particle posi- 
tions and velocities according to the quadrupole formula 



(47) 



As this disk is nearly Newtonian, higher multipoles are 
expected contribute negligibly to h + . The good agree- 
ment over several oscillations in Fig. || tests many facets 
of the mean-field particle simulation scheme, including 
the field and particle integrators, the particle and grid 
sampling, and the source binning algorithm. Wave forms 
extracted at different radii agree at the 15% level at this 
grid resolution. Our experiments indicate that most of 
this difference, and the difference from the quadrupole 
calculation, arises from instantaneously propagated er- 
rors in solutions of the constraint equations. To re- 
move near-zone and gauge effects, the wave extraction 
method requires delicate cancellations between different 
multipole moments. Counter-streaming in the particle 
source can cause inconsistencies between the hyperboli- 
cally and elliptically computed parts of the gravitational 
field, making this cancellation inaccurate. 



C. Cold equilibrium relativistic disk 

Here we consider equilibriu m dis ks constructed accord- 
ing to the prescription of Sec. Ill C . The question we wish 



to answer is whether this disk is unstable to rings as in 
the Newtonian theory when the source and gravitational 
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FIG. 4. Snapshots of the particle positions for the evolu- 
tion of a cold, equilibrium relativistic disk. The initial disk 
radius is Ro/M — 8.0. The rapid growth of concentric rings 
is apparent. 



FIG. 5. Snapshots of the particle positions for the evolution 
of a hot, near-equilbrium relativistic disk. The initial disk 
radius is Ro/M — 10.0 and the velocity dispersion parameter 
Q/flo = 0.85. Following collapse, the disk settles down to a 
new equilbrium state. 



field are evolved consistently. In Fig. [| we show a rela- 
tivistic disk with Rq/M — 8. Before the evolution pro- 
gresses very far, ring formation begins and by t = 30M 
is completely dominant. Less relativistic configurations 
behave similarly. This calculation was carried out with 
300 radial by 16 angular zones and 12000 particles, but 
runs with as many as 48000 particles did not not slow 
down the ring formation. 



D. Hot relativistic disk 



As discussed before, ring formation is suppressed by 
the addition of sufficient random particle motion. Un- 
fortunately, only small values of velocity dispersion can 
be used without dissipating the outer region of the disk. 
Here we consider the evolution of a relativistic disk 
with Rq/M = 10, and velocity dispersion parameter 
f2/f2o = 0.85 and £ = 1.0. Because of its compactness, 
this near-equilibrium relativistic Kalnajs disk is unsta- 
ble to collapse. As shown in Fig. [| the disk initially 
collapses and then oscillates about a new equilibrium of 
about Rq/M = 6.0. At late times it virializcs to a static 
equilibrium state. 

In Fig. H the gravitational wave form extracted at 
r/M = 12.0 is compared with the quadrupole formula 
result. As discussed above, we have evidence that the 
stochasticity of the particle source is responsible for the 
short-time scale discrepancies between the wave forms. 
This calculation used a 200 radial by 16 angular zones 



FIG. 6. Wave form from the hot relativistic disk shown in 
Fig. 5. Wave amplitudes are labeled as in Fig. 3. 
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FIG. 7. Snapshots of the particle positions for the collapse 
of a cold relativistic disk. Initially the radius is Ro/M = 1.5 
and the particles are all at rest. The apparent horizon (dashed 
line) first appears at t ~ 4.0M. 

and 12000 particles. Improving the accuracy requires 
higher grid resolution and more particles; unfortunately 
the noise in the wave form is a y/N effect so increasing 
the number of particles quickly becomes computationally 
prohibitive. 

E. Cold disk collapse 

Finally, we consider the collapse of a cold homogeneous 
disk, the disk analogue of Oppenheimer-Snyder collapse 
in spherical symmetry. Initially, all particle velocities are 
equal to zero (£ = 0). We consider a very relativistic case 
with Ro/M = 1.5. (Recall that in isotropic coordinates 
a Schwarzschild black hole has a radius Ro/M = 0.5.) 
With such a compact disk, the collapse is quick enough 
that an apparent horizon appears before the ring instabil- 
ity becomes significant. In Fig. |?] we show results from an 
evolution carried out with 300 radial by 16 angular zones 
and 24000 particles. The apparent horizon appears at a 
time of about 4.0M and ring formation is just discernible 
at the disk center at this time. 

In order to prolong the numerical evolution after the 
black hole forms and counteract the effect of "throat 
stretching" that occurs in our maximal (and other sin- 
gularity avoiding) time-slicing conditions, we have imple- 
mented a moving mesh algorithm that moves the inner 
radial grid-zone to track the growth of the conformal fac- 
tor i/j. With this method, we are able to evolve the black 
hole and preserve constancy of out quasi-local mass mea- 
sures (the Brill and ADM masses) to a few percent for 
around 15M after hole formation. Evolution beyond this 
time is prevented by the numerical difficulties in inte- 



FIG. 8. Gravitational wave form from the disk collapse 
shown in Fig. 7. The quadrupole wave form extracted at 
r/M = 6.0 is shown as a function of retarded time. 

grating the particle geodesies on the extremely stretched 
mesh. 

This time is sufficient, however, to see the pulse of 
radiation from collapse and, perhaps, the first half wave- 
length of a quasi-normal mode oscillation. In Fig. |§| we 
show an estimate of the asymptotic quadrupole wave 
form extracted at r/M = 6.0. After an initial peak repre- 
senting radiation in the initial data we see an oscillating 
signal with total wavelength of approximately 16Af , com- 
parable with that of the most slowly damped 1 = 2 quasi- 
normal mode oscillation which has A = 16. 8M. The en- 
ergy radiated is less than 0.1% in the quadrupole mode 
and we estimate that it is less than 0.01% in higher l- 
modes. 

We monitor the area of the apparent horizon and its 
polar and equatorial circumferences (see Ref. |6) for def- 
initions). As shown in Fig. ^, the normalized apparent 
horizon area ^4/167rM 2 has a value of 0.95 when it first 
forms, and asymptotes towards the Schwarzschild value 
of 1.0 before the calculation terminates. The black hole 
is initially oblate when the apparent horizon forms. The 
normalized proper circumferences oscillate around their 
Schwarzschild values C m /AmM — C p /AttM — 1.0 as the 
black hole radiates off its initial asphericity. This oscil- 
lation period again is comparable with that of the most 
slowly damped 1 — 1 quasi-normal mode. 
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FIG. 9. Horizon diagnostics for the black hole formation 
shown in Fig. 7. The proper equatorial circumference, po- 
lar circumference, and area of the apparent horizon normal- 
ized to their Schwarzschild values, C^/AirM, Cp/AuM and, 
A/l&itM 2 are plotted as functions of time in units of M. 



V. CONCLUSIONS 

We have presented here preliminary results from a new 
relativistic mean-field, particle simulation code that can 
evolve disks of collisionless matter and compute the emit- 
ted gravitational radiation. The disk analogue to the 
Oppenheimer-Snyder solution, cold homogeneous disk 
collapse, is an interesting benchmark for numerical rel- 
ativity codes as it is one of the simplest scenarios that 
encounters the two most computationally challenging fea- 
tures of relativistic collapse: black hole formation and 
evolution, and gravitational wave production and prop- 
agation. The hydrodynamical formulation of this prob- 
lem, provided in Sec. II C, 



should permit implementa- 
tion of disk collapse in codes without collisionless matter 
sources. 

Failure to complete the computation of the full radi- 
ation data in the case of disk collapse to a black hole 
reflects the fundamental problem in numerical relativ- 
ity. Specifically, no general algorithm is currently known 
that can integrate Einstein's equations for long enough 
after a black hole forms to compute the full gravitational 
wave form. It may be possible that this problem can be 
solved with suitable apparent horizon boundary condi- 
tions, "cutting out" the black hole (see e.g. Ref [0). 
We expect disk collapse to be a useful proving ground 
for developing apparent-horizon boundary conditions and 
other black hole evolution techniques in an axisymmetric 
context. 
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TABLE I. Disk Evolution Cases 



Case 


Ro/M 


n/n 




fate 


Oscillating cold Newtonian (analytic source) 


30.0 


1.0 


0.9 


oscillates 


Oscillating Kalnajs 


20.0 


0.8 


0.9 


oscillates 


Cold relativistic equilibrium 


8.0 


1.0 


0.99 


rings 


Hot relativistic near-equilibrium 


10.0 


0.85 


1.0 


collapses and virializes 


Cold relativistic 


1.5 


1.0 


0.0 


collapses to black hole 
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